Tunneling effects in a one-dimensional quantum walk 



O 

(N 

Oh- 
< 

in 

(N 



a,: 

-I— > 

c ■ 

ctf ■ 
=5 

cr 



> 

in 

m 

o 
o 



X 



Mostafa Annabestani/'B Seyed Javad Akhtarshenas, 2 ' 3 '0 and Mohamad Reza Abolhassani 1 

'Department of Physics, Basic Sciences Faculty, Tarbiat Modares University, Tehran, Iran 
2 Department of Physics, University of Isfahan, Isfahan, Iran 
'Quantum Optics Group, University of Isfahan, Isfahan, Iran 
(Dated: April 27, 2010) 

In this article we investigate the effects of shifting position decoherence, arisen from the tunneling 
effect in the experimental realization of the quantum walk, on the one-dimensional discreet time 
quantum walk. We show that in the regime of this type of noise the quantum behavior of the walker 
does not fade, in contrary to the coin decoherence for which the walker undergos the quantum-to- 
classical transition even for weak noise. Particularly, we show that the quadratic dependency of the 
variance on the time and also the coin-position entanglement, i.e. two important quantum aspects of 
the coherent quantum walk, are preserved in the presence of tunneling decoherence. Furthermore, 
we present an explicit expression for the probability distribution of decoherent one-dimensional 
quantum walk in terms of the corresponding coherent probabilities, and show that this type of 
decoherence smooths the probability distribution. 

PACS numbers: 03.67.-a, 05.4a.Fb, 03.65.Xp 



I. INTRODUCTION 

The quantum walk (QW) is the quantum version of the 
classical random walk (CRW). The superposition prop- 
erty of quantum systems allows us to define quantum 
walker which can walking in all possible paths, simulta- 
neously, and the interference of these paths makes some 
important differences between the QW and the CRW. 
Notable differences between the QW and the CRW are 
the quadratic dependency of the variance on the number 
of steps and the complex oscillatory probability distri- 
bution in the QW instead of the linear variance depen- 
dency and the binomial probability distribution in the 
CRW. These differences have been used to present sev- 
eral quantum algorithms in order to solve some specific 
problems [lHg] with performances better than the best 
known classical versions. 

Two types of QW have been introduced: discrete time 
Q and continuous time (sj . Although the mathematical 
structure is completely different in these two types, but 
both have the same results in our problems. The relation 
between these two types of QW was an open problem for 
several years, but this relation has been found recently 
9, 10]. A lot of works have been devoted to the study 
of QW, both from theoretical as well as experimental 
view points. Many aspects of QW have been studied 
theoretically, for example some researches focus on the 
network over which the walks takes place. In this area, 
the QW on a line has been well studied 

EMU, 

but other 

topolog ies such as cycles [13, EU, two-dimensional lattices 
[l~9l - l22| ]. or n-dimensional hypercubes [H, [24| have also 
been investigated. 

The study of quantum entanglement and its connec- 



tion with the QW is another important research area 
which has recently attracted much attention. The ef- 
fect of entanglement of the coin subspace on the QW 
[3 HH HEI , entanglement between the coin and the posi- 
tion subspaces I20l |26|. 12711 and the QW as the entangle- 
ment generators [2g, l29j are examples of these studies. 

Beside of all theoretical studies, the experimental im- 
plementation and realization of the QW are also re- 
ceived much attention [3fj| - l35| . In the experimental im- 
plementation, the environment effects will be so impor- 
tant because, in practice, the preparation of pure quan- 
tum states without interaction with the environment is 
impossible, and the environment can disturb the quan- 
tum states and fades the quantum properties. Therefore, 
it is very important to formulate and quantify the in- 
fluence of decoherence on the QW and, indeed, several 
valuable works have been made in this subject [32l |36| - 
44j |. Moreover, the study of decoherent QW can be con- 
sidered as a link between the theoretical study and the 
experimental implementation. 

It is known that the quantum properties of QW, such 
as the quadratic dependency of variance on the time, 
make the QW more powerful than the CRW, but these 
properties are very fragile when they are exposed to de- 
coherence or noise. Brun et al. [3!| have shown that the 
quadratic term of variance vanish, in long time limit, even 
for weak decoherence of the coin subspace. On the other 
hand, the linear term remains greater than the classical 
one for reasonable strength of noise. It means that if the 
noise on the system be less than some certain value, the 
QW will spread faster than the CRW. Since noise on the 
quantum systems is unavoidable and is, indeed, a mecha- 
nism for emergence of the classical behavior in quantum 
systems, some authors have tried to use such classical be- 
havior emerged by the noise. For example, Kendon et al. 
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37| have shown that a weak noise, both on the coin and 
the position subspaces, can be useful in quantum algo- 
rithms. They have explained that a general weak noise 
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can smooths and flats the probability distribution. Of 
course, we should note that this smoothness occurs in 
the cost of losing the speed of spreading (variance). On 
the other hand, other route to classic behavior of quan- 
tum walk has been investigated and compared to classical 
behavior emerging by decoherence. Brun et al. [39j ] have 
shown that multiple coins version of QW retains "quan- 
tum" quadratic growth of the variance except in the limit 
of a new coin for every step. 

In this paper we use our analytical expression for vari- 
ance, recently presented in [44j, and show that the shift- 
ing position decoherence does not change the quadratic 
behavior of the variance in the 1DQW. We show also 
that although this quadratic term depends on the initial 
state, but it can not be removed even for initial mixed 
states. Furthermore, we show that the effect of this type 
of decoherence is to smooth the probability distribution 
and we present a closed formula for the probability distri- 
bution in terms of the probabilities of the corresponding 
coherent QW. 

This work is organized as follows. Section II gives a 
brief review on the 1DQW and decoherence. Our calcu- 
lations and results for tunneling effects have been put in 
section III. In this section we obtain an expression for the 
variance of 1DQW in the presence of tunneling noise and 
show that the time quadratic dependency of the variance 
exists for all initial states. The entanglement properties 
between the coin and the position degrees of freedom 
have been also studied in section III, and it is shown 
that in this case, unlike to the case that only the coin 
is subjected to decoherence, the entanglement converges 
to a significant value. Moreover, a closed formula for 
the probability distribution has been presented in this 
section. We summarize our results and present our con- 
clusions in section IV. 



II. BACKGROUND 

QW is the quantum version of the CRW where instead 
of coin flipping, we use the coin operator to make super- 
position on the coin space, and instead of the walking 
we use the translation operator to move quantum parti- 
cle according to the coin's degrees of freedom. In one- 
dimensional quantum walk we have two degrees of free- 
dom in the coin space H c , spanned by {\L), \ R)}, and infi- 
nite degrees of freedom in the position space H v , spanned 
by {\i) i = — oo, ■ • • ,00}. The walker Hilbert space Hyy 
is defined as the tensor product of the coin space H c and 
the position space H p , i.e. H\y = H s <E> H c . Each step of 
QW is constructed by the unitary operator U c , making 
a superposition on the coin space, followed by the trans- 
lation operator 5, moving the particle according to the 
coin state, i.e 

U W =S(I®U C ), (1) 



where 

S = N + 1) <a;| ® \R) (R\ +\x-l) (x\ ® \L) (L\. (2) 

X 

Therefore, the unitary evolution of the quantum walker 
is as follows 

|*(t + l))=E/ w |*(t)) \*(t)) = U v> t \*(0)). (3) 

This is, in fact, the coherent evolution of the system but, 
in practice, it is not possible to isolate the system from 
the environment and, in general, the environment affects 
the coherent evolution of the system. Generally, deco- 
herence is used to estimate deviation from the ideal case 
in which the effects of environment are neglectable. One 
important approach to investigating decoherence is the 
so-called Kraus representation (45[. Let us define He as 
the Hilbert space of the environment, spanned by {|e n ), 
n = • • • m} where m is the dimension of the environ- 
ment's Hilbert space. Now, in order to study the time 
evolution of the system we should consider the time evo- 
lution of the whole system, i.e. system+environment, 
defined on the Hilbert space H = He ® %, and obtain 
the state of the system by tracing out over the environ- 
ment's degrees of freedom, i.e. 

p sys = Tr Env (UpU*). (4) 

Here U acts both on the system and the environment 
Hilbert spaces. Without loss of generality we assume that 
the state of the whole system is p = po <g> |era;o)(era>o|. 
So we can write Eq. Q as 

m m 

Psys = ^2 (e„| U \env ) p (env \ |e„) = ^ E n p El, 

n=0 n=0 

where E n = (e n \ U \envo) , n = 0, 1, • • • ,m are the so- 
called Kraus operators. These operators satisfy the fol- 
lowing completeness relation 

J2 E * E n = I' (6) 

71 

By definition of the Kraus operators, one step of the walk- 
ing can be written as follows 

m 

p{t + l) = Y J E n p(t)El (7) 

n=0 

For t steps, we can write 

m m m 

p(t) = E - E E E nt ...E n2 E ni p(0)ElEl 2 ...El. 

n t =0 n2=0 ni — 

(8) 

It is worth noting that Eq. ([8]) is general and the Kraus 
operators E n include the whole information of all types 
of evolution. It follows, therefore, that the coin operator, 
the translation operator and the environment effects all 



3 



are embedded in E n and we did not assume any restric- 
tion yet. 

To make any progress, we should therefore find the 
Kraus operators for our system defined in Eq. ([5]), and 
use Eq. ([8]| in order to obtain the final state p(t). Ev- 
idently, the E n are operators that act on the system 
(coin+position) Hilbert space, and therefore we can write 
the general form of E n as follows 

E n = EE a ^«^')^®l l )0'l ( 9 ) 

x I i,j 

where x, I = — oo, • • • , oo and i, j — {L, R}. 

Recently, we have shown that a reasonable suggestion 
on the coefficient a^- n ' of Eq. © enables us to derive 
useful analytical expression for the first and second mo- 
ments of position [44f j . In the remaining of this section, 
we briefly introduce the method. To begin with, we use 
the Fourier transformation 

7T 

-e-^|fc), (10) 

— 7T 

and write Eq. (|9j) in the /c-space as 

& Y.Y.<:.., I f ^ e -«*e- fa (*-*') |fc> (fc'| ® \i) 0'|. 

x,l i,j 

(11) 

If we assume that the coefficients a^J . . are not de- 
pendent on the coordinate x, then one can use the or- 

thonormalization relation e~ M ( fe ~ fc ') = 2-k8 (k - k') 

X 

and caries out the summation on x. This is, of course, a 
reasonable assumption for a large family of QW, includ- 
ing the QW on the homogeneous position space. With 
this assumption, we have 

En = / l fc ) ® C n (*)> (12) 

J 27T 

where 

Cn(k) = J2I2 a i% e ~ Uk \ i )(i\- ( 13 ) 
Now, if we write the general form of po in the fc-space as 
Po= J J -^\k)(k'\ ®|Vo) (Vol, (14) 

we get, from Eq. (J7J, the following form for the first step 
of walking 

f f rjhfjh 1 , 

P = ^\k){k'\®Y. c n (k)\^){MC n (k') 

J J n 

= J J ^r\k){k'\®L k , kl |Vo) (Vol, (15) 



where the superoperator Ck,k' is introduced by 

C k , k ,6 = Y J Cn{k)dCi{k'). (16) 

n 

Therefore, after t steps walking, we find 

p(t) = J J ^rlfc>(fc'l®4, fc ' IVo) (Vol, (i7) 

as the state of the walker and 

p(M) = jJ^(x\h)(h'ix)Tt (18) 

= //^«-('"-'>T r (4,«). 

as the probability of finding the walker in the position 
x. Note that the completeness relation on the Kraus 
operators, given in Eq. ©, now implies that the coin 
operators C n (k) satisfy the same relation, i.e. 

J2C'Uk)C n (k)=I, (19) 

n 

where can be used to prove another important property 
of Ck,k, i-G. the trace preserving condition 

Tr (Cl k 6) = Tr (6) , (20) 

for arbitrary operator O. 

By definition, the mth moment of the probability dis- 
tribution p (x, t) is as follows 

(x m )=J2^PM, (21) 

X 

where can be used , together with Eq. (TT51) . and write 
the first and second moments as 

W = hff ^'f^O^Ai) (22) 

Now, by using the relations 

^-Tr(£ kk ,p)=Tr(g kk ,p) (23) 
ak 

-^Tr [Ckk'P) = Tr (G f k >kp) , 

where 

g Kk ,d = Y,^^ocUn (24) 

n 

one can carries out the integrations of Eq. (|22|) and, after 
some calculations, get the following relations for the first 
and second moments 
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(x) t = i 
(A = 



dk 



(25) 



* dk 



m — 1 
t m—1 



2tt 



m—1 m' — l 
m — m —1 / ^-jf ^m 



+ 



0JC -1 liM<*>l)} 



Here we have defined, for simplicity, Lk = Ck,k, Qk 
Gk,k and 



Jk 



dQlk' 



dk 



k'=k 



dk 



dk' 



(26) 



k' = k 




10 20 30 



FIG. 1: (Color online) Coherency of system versus t for var- 
ious p. Note that for p — (coherent case), the coherency 
does not decay. 



III. CALCULATIONS AND RESULTS 

In this section, by using the analytic expressions for 
the first and second moments given in Eq. (|25[) . we in- 
vestigate the variance of the 1DQW in the presence of 
tunneling noise, i.e. a type of decoherence which is ap- 
plied on the position subspace. 



A. Variance of decoherent 1DQW 

Consider a one-dimensional QW such that after each 
step of walking, the walker can move to the nearest neigh- 
bors with probability p. This phenomena occurs in the 
experimental setup because of tunneling effect. The ef- 
fect of such noise on the probability distribution has in- 
vestigated, numerically, by Diir et al [32| ■ In this section 
we calculate, analytically, the variance of 1DQW in the 
presence of this kind of noise and investigate some im- 
portant features of it. 

Let us assume that, after each step in 1DQW, the 
walker tunnels to the left or to the right site with prob- 
ability p, so it is clear that the walker does not move 
with probability (1 — p). Therefore, unlike the coherent 
QW, all sites may be occupied. In the presence of such 
tunneling, one step of the walking can be written as 



P (t+i) = (i-p)u wP (t)ui 

+ | (s + U wP (t) Utsl + S-U wP (t) uts: 



(27) 



I r and U,„ is the coherent 



where S± = ^\x±V) { 

X 

walking operator, given by Eq. (fTJ). For p = 0, the 
system is exactly the same as the coherent QW and the 
evolution of the system is unitary, so the final state of 
the system remains pure. On the other hand, for p =/= 
the evolution is nonunitary and decoherenc is happened. 

Generally decoherence can be considered as a sym- 
bol to express the reduction of purity of the system 



state and, therefore, one can use the linear entropy 
S(p(t j) = 1 — Tr[p 2 (t)} as a measure of decoherence. The 
linear entropy has the limiting values and 1 — 1/N, re- 
spectively, for pure and maximally mixed states, where 
N is the dimension of the space that the density matrix 
p(t) is supported on. In the same manner, one can use 
the purity (or coherency) 



CW=Tr[p 2 (t)] 



(28) 



where has the limiting values 1 and 1/N for pure and 
maximally mixed states, respectively. Therefore, the pu- 
rity less than 1 is a signature of decoherent evolution. In 
Figll] we have plotted the coherency of the system as a 
function of time in the presence of the tunneling noise of 
Eq. (|27p. for different values of the probability p. As we 
expect, for p = the evolution of QW takes place coher- 
ently, while by increasing p the system loses its coherency 
very fast. 

Now by comparing Eq. (|27[) with Eq. (J7J , we can write 
the Kraus operators as 



(29) 
(30) 




(31) 



All the above Kraus operators have the form given by Eq. 
(|9]), therefore, by using Eq. (fF3|) we find the Ci matrices 
as 



_ ,'l - p f e lk e lk 



(32) 



l -l 



r _ Vp ( 1 1 

°3 — — I £ 2ik _ e 2ik 
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where satisfy, clearly, Eq. (|19[) . 

Now by using the afhne map [45| , we can represent ar- 
bitrary two-by-two matrix O by a four-dimensional col- 
umn vector as 



and 



O = ^ r ^ a i 
i=0 



(33) 



where r, = Tr(aiO)/2. Now, in order to find the super- 
operators CkiGkiGf. and Jhi we use the affine map and 
represent the action of them on an arbitrary two-by-two 
matrix O as follows 



C k O 



Q k O 



( 1 

sin(2fc) cos(2fc) 

-cos(2fc) sin(2fe) 

\0 1 





(r \ 




ri 






) 


\r 3 J 







—i 



-i 

cos (2k) 

sin (2k) 





\ 

- sin (2k) 

cos (2k) 




/1 + 



J k O 



1 -hp 








qsin(2fc) q cos (2k) 

— gcos(2fc) gsin(2fc) 

1 +p 



fro\ 
r\ 

\r 3 J 

\ 



(34) 



(35) 



fro\ 

ri 

\r 3 J 



(36) 

and Qj. = C/j*, which is obtained from the Hermiticity of the 
Pauli matrices. 

Now, according to Eq. ([251) . in order to calculate mo- 
ments and variance, we need to find the mth power of 
Ck- But, here, the method of non-trivial submatrix 
which is introduced in [38], is not applicable because 
det (/ — M) = 0, i.e. I—M is not invertible. Fortunately, 
the eigenvalues and eigenvectors of Ck are enough sim- 
ple to calculate C™ directly. Straightforward calculations 
show that the eigenvalues and the corresponding eigen- 
vectors of matrix Ck, given in Eq. (|3~4")1 , are as follows, 
respectively 



Ai — A2 



1 



A, 



ol (e+7r) 



A 4 



-i(0+7r) 



(37) 



ei 



1 



l e 2) = -± 



C3 



sin (k) \ 
cos (k) 
sin (k) 
cos (ft) / 
-2 (cos 2 (ft) + l) \ 

sin (2k) 
2 sin 2 (k) 

sin (2fc) / 

/ 

sin(2fc) e 2i " 

2sin 2 (fc)e 2lS 
V - sin (2k) e i6 



(38) 



1 

N 3 



where cos (9) = cos 2 (k) and iVj are the normalization 
factors. By the help of eigenvectors and the spectral 
decomposition, we are able to find any power of matrix 
Ck, appeared in Eq. ([251) . 

In order to calculate the first moment (x) , we need 



T = ^2cT- 1 = ^2{\ei){e 1 \ + \e 2 ){e 2 \ 



m — 1 

+ei(e+x)(ro _ 13 |es) (ea| + e _ l(9+ , )(m _ 1) |e4) ^ j 

t(|ei)(e 1 l + |e 2 )(e 2 |) 

'l_ e «(«+»)\ /i_ e -i*(9+-)\ 
+ I : . s 1 |e a > (e 3 | + - ■ _.. |e 4 ) (e 4 | 



1 + e»« 



1+e 



= * (|ei) (e a I + |ea> H) + ™(j^e ^ H 



(39) 



Here, in the last equality, we used the fact that the forth 
eigenvalue and eigenvector are complex conjugate of the 
third one, respectively (see Eqs. (|37|) and (|38])1. Also we 
neglect the oscillatory term e lt ^ e+ ' R \ because according 
to the phase stationary theorem the oscillatory term can 
be neglected in long time limit. 

By putting Eq. ([38]) into Eq. (|39|) and some simplifi- 
cation, we find the following matrix form for T 



m — 1 



1 

A 



/ 2t [ 



< os- (k) + l] 






2£cos 2 (k) + 1 
(i - 1) sin (2Jfe) 
2 (t - 1) cos 2 (ft) 4 





t sin (2ft) 
cos 2 (Ar) + t sin 2 (k) 
(t - 1) sin (2ft) 



\ 

2icos 2 (k) - 1 

t sin (2ft) 
2icos 2 (ft) + 1 ) 



(40) 



where A = 2 [cos 2 (ft) + l] . We can, therefore, write the 



first moment (x), given in Eq. (|25|) . as follows 



W« = 



2tt 



dfc ( -2i ) r 



ri 

^2 
f"3 



(41) 



(2-V2)t + i_ 



ri 



V 



'"3 • 



6 



Note, however, that in calculating this moment from Eq. 
(f2"5| only the first row of Qk gives the nonzero contri- 
bution. Clearly, the above expression for (x)t does not 
contain p, which implies that this type of noise can not 
change the first moment. 

Now, in order to calculate (x 2 ), we first calculate the 
last term of Eq. ([25)1 as follows 



On the other hand, the first term of Eq. (|25[) can be 
written as follows 



1 

2^ 



rffc^Tr^^r'lV'o)^!)} 



dk ( i + p o o o ) r 




t(l+p). 



(42) 



/: 



m— 1 m'—l 



^ E E Tr {sic?-™'- 1 (QkCf- 1 |Vto)<^|) + GkCT- 1 '- 1 (gIct'- 1 \rhHM) } (43) 

-/ *EE(° i00 ) £ r m " 1 (&-^r" 1 

•' _7r m=lm'=l 




where we have used the fact that Q\ — Qt- From Eq. where V is defined by 
([35]) the exact form of £/& - Q\ is t 






-2i 


















!) 

















-2i 













Also from the trace preserving property of we have 



So, easily, we get 





/1/2 












V ^ 




/- 



(45) 



(46) 



Putting this into Eq. (|4*31 . we get 

rn dk 



f *(oioo)r' 



(47) 



ei) (ei + e 2 ) (e 2 



m=l m'=l m—lm'=l 

(44) +e i(0+^)(m-m'-l) ^ ^ + e -<(9-Hr)(tn-m'-l) ^| J 

= ~(*-l){|ei>v3i| + |e 2 )<e 2 |} (48) 

-2 !^ I ) lea) (e 3 | 1. (49) 



v x 1 + e (l + e' ie ) 7 , 

From Eq. (|47)l it is clear that we need only the term T' 2 4 , 
which is 



r' 



t 2 (cos 4 (k) + cos 2 (fc)) + 1 t 
2 (cos 2 (jfej + l) 2 2' 



(50) 



Using this into Eq. (|47|) . caring the integration and 
putting everything together, we find the following form 
for the second moment 



(x 2 ) = t(l+p)+ 1- 



1 

71 



t 2 -t 



3-y/2 



1 



V2 



2 3\/2 
t 2 



(51) 



It is worth to note that the second moment does not 
depend on the initial state and, surprisingly, it involves 
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the quadratic term t 2 . This means that the decoherent 
1DQW preserves the quantum property (ballistic behav- 
ior of variance) whenever the position space is subjected 
to decoherence (tunneling noise), in contrary to the case 
that the coin is subjected to noise in which for long time 
walking the system becomes classic and its variance be- 
comes linear, even for weak strength of noise [381 ] . 

Let us look at the variance more precisely. From Eqs. 
dSJ and ([51]) we have 

V = (x 2 ) - (x) 2 = At 2 + Bt + C, (52) 

where 

A = a - Aa 2 (r 3 + r x ) 2 

B = 2V2a(rl-rf) + p 

C = -(r 3 -r 1 ) 2 /2 + 3^2/8, (53) 

in which a = 1 — l/v2- In the long time limit and for 
A 0, we can neglect the B and C coefficients. The 
question that may arise here is that: is it possible to find 
initial states for which A = and B ^= ? or, in other 
words: is it possible to find initial states such that the 
variance is classic?. It is not difficult to show that it is 
impossible to have A — even for mixed initial states. 
Let us consider the following general form for the initial 
state po 

E( r +r 3 ri-ir 2 \ 
(Tin = i • , (54) 

\ri+ir 2 r -r 3 I 1 y ' 

i=0 v 7 

where ro = 1/2, because of the normalization, and the 
positivity of po implies that 

r\+r 2 +r 2 < ± r\ + r 2 < J. (55) 

This means that \rx\, \r 3 \ < 1/2 or 1r\r 3 < 1/2, so 

(ri+r 3 ) 2 =r 2 +r 2 + 2 ri r 3 < |. (56) 

On the other hand, the coefficient A in Eq. ([53")) vanishes 
only for 

(ri + r 3 ) 2 = « 0.85, (57) 
4a 

which is, obviously, in conflict with the condition given 
by Eq. (|56[) . Consequently, the tunneling noise which is 
a noise on the position subspace of QW, not only dose 
not change the quadratic behavior of variance (p does 
not appear in the coefficient A in Eq. (|53p). but also 
increases it a bit (as a linear term in the coefficient B). 
Furthermore, although the quadratic term depends on 
the initial state, but it can never be vanished by tuning 
the initial state and it ranges between the minimum a — 
2>a 2 and the maximum a. Now, in order to find the initial 
state for which the variance is maximum, let us consider 
the following generic form for the pure initial state 

|V»o> =cos6\R) + e < *sin0|£) ) (58) 



800 
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FIG. 2: (Color online) Variance versus £ for different initial 
states of Eq. (158[) . The continues curves are the results of the- 
oretical formula (|52[) and the corresponding discrete symbols 
are the direct numerical calculations. 

where < 8 < tt and < <j> < 27r. For this initial state, 
the elements of in affine map representation defined in 
Eq. (|33f are as follows 

ro = I r\ = i cos sin 2$ , , 

r 2 = \ sin 4> sin 26 r 3 — | cos 29. 

Using these r\ , r 3 in the definition of A, given in Eq. (|53|) , 
and maximizing it, we find that for all 9, <fi, satisfying 
(f> = — cos" 1 (cot 29), the variance is maximum. In Fig. 
[21 we have plotted variance versus time for various initial 
states and have compared the results with the numerical 
simulation. The figure shows that there are excellent 
matching between the theoretical result of Eq. ([52[) and 
the direct numerical calculations. An interesting point 
which we would like emphasize here is that, although 
our expression have derived for long time limit, but very 
fast decay of oscillatory term implies that we have very 
good matching even after a few steps (see Fig. ©). 

B. Entanglement in decoherent QW 

In the previous subsection we showed that the tunnel- 
ing noise preserves the quantum property of variance. In 
this subsection we turn our attention to another quan- 
tum property of the system, namely, the entanglement 
between the coin and the position degrees of freedom, in 
the presence of tunneling noise. We also compare the 
result by the case that only the coin is subjected to de- 
coherence. The model that we choose for the coin deco- 
herence is defined by 

p(t + l) = (l-p)U wP (t)Ul (60) 
+ P [PrU wP (t) UlP R + P L U wP (t) UlP L ] , 

where Pr = I X ®\R) (R\ and Pl — I x ® \L) (L\ are projec- 
tions on the coin subspace. Equation (pO)) describes the 
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situation in which, after each step of walking a measure- 
ment can perform on the coin subspace of the walker with 
probability p. The probability distribution and the vari- 
ance of this type of decoherence have been investigated 
inRef. 

As we mentioned before, in the presence of noise, the 
state of the system becomes mixed and, therefore, we 
use the negativity of the partial transpose of the den- 
sity matrix as a suitable measure of entanglement which 
is, of course, a computable measure for mixed bipartite 
systems [46| • Negativity is based on the Peres- Horodecki 
criterion for separability [U |48[ and measures the degree 
to which the partial transpose of p fails to be positive, 
i.e. the absolute value of the sum of negative eigenvalues 
of the partial transpose of the density matrix 



N = 



(61) 



where A^ and A^ are eigenvalues and negative eigenval- 
ues of the matrix p'(t). Here p'(t) is the partial transpose 
of p(t), which is obtained from the density matrix p(t) 
by taking the transpose of p(t) with respect to one of the 
subsystems, say the first subsystem, i.e. 



Pxc,yb Pxb,yci 



(62) 



where x, y and 6, c denote the position and the coin state 
indices, respectively. 

In Fig. © , we have plotted negativity of the system 
in terms of the noise strength p, both for the cases that 
the position is subjected to the noise given in Eq. (|27|) 
and the coin is subjected to the noise given in Eq. (|60| . 
Comparing two negativities, it is clear from the figure 
that the system with coin only decoherence is more in- 
fluenced by the noise and losses its negativity very fast. 
Moreover, though the negativity of the system with coin 
only decoherence decays to zero value as p grows up, but 
situation is very different for the system with position 
only decoherence and, indeed, for such system the neg- 
ativity converge to a significant value, never less than 
« 0.7, even for maximum rate of noise. This means that 
the system with tunneling noise never loses all quantum 
behavior, but the system with coin decoherence will do. 
In Fig. (@}, we have plotted negativity of the system as 
a function of time, in the presence of tunneling noise for 
various p. The figure, clearly, shows the convergence be- 
havior of the negativity for this type of noise. Note that 
the flat behavior of the negativity in coherent case, i.e. 
p = 0, changes to descending (ascending) behavior for 
p < 0.5 {p > 0.5) (see Fig. gj). 



C. Probability distribution of decoherent QW 

In this subsection we derive an analytical expression 
for the probability distribution of decoherent quantum 
walk in terms of the corresponding coherent probability 




- Position only decoherence 

- Coin only decoherence 



FIG. 3: (Color online) Negativity of the system with respect 
to t in the presence of the position only decoherence Eq. (|27[) 
(thick curve) and the coin only decoherence Eq. (|60[) (thin 
curve) . 
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FIG. 4: (Color online) Negativity of the system versus t, in the 
presence of tunneling noise. The curves from top to bottom 
show the negativity when the noise strength p goes from to 
1 with 0.1 step length. The upper curve (thick curve) is the 
negativity of the coherent case, i.e. p = 



distributions. To begin with, we first rewrite Eq. (|18l) in 
the coherent case p = as 



P (x,t) 



dkdk' 

47T 2 



e-^ fe '- fe )Tr(w£ fe ,p ). (63) 



Note that we use Po (x, t) and Wfc,fc' instead of P (x, t) 
and £k,k' , respectively, in order to denote the correspond- 
ing quantities for coherent QW. It is clear from Eqs. (fT6|) 
and (El that 



Wk,vPo = U{k) polJl (A/), 



(64) 



where U (k) is the Fourier form of the unitary transfor- 
mation of the Hadamard walk in the fc-space, i.e. 



U(k) 



V2 ( 



(65) 
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Accordingly, all given in Eq. (|32l can be rewritten as 
follows 



C 2 = 



(66) 



C 3 



%e- %k U{k) 
%e tk U(k). 



We can , therefore, write 
£ k ,k>Po = ^C n (k) Po Ct(k') 



[l_ p+ P( e <fc'-fc) + e -<fc'-fc) 



(67) 



Putting this into Eq. (fT8[) and using Eq. (l63l) . we find 
the following form for the probability distribution of de- 
coherent QW 



n— m— 



'co-EE UJur (68) 

(l-p)'""(|)"Po(!» + 2m-n,«), 



where 



is the binomial coefficient. 



Equation (1551) shows that, after i steps of decoherent 
walking, the probability of finding the walker at site x 
can be written as a linear combination of the probabili- 
ties of the corresponding coherent walking, but at some 
other sites ranging from x — t to x + t. The equation 
can be also used in order to prove the smooth property 
of the distribution and all sites occupation property. A 
comparison of this formula and the numerical simula- 
tion shows that there are an excellent consistent between 
them. Figure ([5]) shows the probability distribution in 
terms of x, in the presence of position decoherence ( tun- 
neling noise ) with different noise strength. It is worth to 
mention that although Pq (x, t) = for odd (even) sites 
after even (odd) steps, but it is evident from Eq. ([55)1 
that the probability P (x, t) of finding the walker in all 
sites x and < p < I is nonzero (see Fig. ((5}). In the 
totally decoherent case, i.e. p = 1, only the terms with 
n = t have nonzero contribution in the summation of Eq. 
|, leads therefore to 



P(x,t) 



)=i 



I 4 

2* 

m=0 



P Q (x + 2m - t, t) (69) 



This equation shows that in the case p = 1, unlike the co- 
herent case p — 0, only even sites are occupied for any t. 
The reason is that if t = 2k, then only for x = 21 we have 
P(x + 2m-t,t) = P (2s, 2k) ^ with s = I + m - k. 
On the other hand, if t = 2k + 1, again only for x = 21, 
we have P (x + 2m - t,t) = P (2s - 1, 2k + 1) ^ 0. This 
means that the walker influenced by full tunneling noise, 
can not be found in odd position after any number of 




FIG. 5: (Color online) Probability distribution of the position 
only decoherence (thick curve) for (a) p = 0.01, (b) p = 0.1 
and (c) p = 0.9. The thin curve is the probability distribution 
of the coherent QW, which is plotted for all sites and allows 
us to compare all site occupations in the position decoherent 
QW with odd (even) site occupations in the coherent case 
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FIG. 6: (Color online) Probability distribution for strong val- 
ues of noise (a) p = 0.99, (b) p = 0.995 and (c) p = 1. The 
thick and smooth curve is, in all cases, the probability dis- 
tributions for p = 0.97 which is approximately critical point 
in which the behavior of the probability distribution changes 
from smoothness to rough distribution by growing p. 



steps. Numerical calculations show that for p ^ 0.97 the 
probability distribution loses it's smoothness and con- 
verges to the case p = 1 in which only half of the sites, 
namely even sites, are occupied. Figure ([6]) shows the 
probability distribution for some large values of p. It is 
evident from Figs. ([5]) and © that for very weak and 
very strong noises the probability distribution is rough, 
but for the other strengths of noise the tunneling effect 
can smooths the probability distribution. In order to 
quantify the smoothness of the probability distribution 
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in terms of p, we use total variation of the probability 
distribution as a measure of roughness of distribution, 
defined by 



Var» (/) = sup 



ns — 1 



8=0 



\f(Xi+l) - f(Xi)\, 



(70) 



where the supremum is taken over the set S = 
{S = {xq, x ns } |, S is a partition of [a, b]}. This def- 
inition is calculated in the interval of all possible sites 
and is a relevant measure for quantifying roughness of 
the probability distribution. In our problem we have 



Var*_ 4 (P) = J2 \P(i + l,t)-P(i,t)\. 



(71) 



Note that in interval [— t, t], the set S = {— t, ...,£} in- 
cludes all sites and we have discreet values for probability, 
so supremum occurs automatically. 

Figure (JTJ) shows the total variation of the probability 
distribution with respect to p. As we expect, for small 
value of p, the total variation is large, which means that 
the probability distribution is not different from the co- 
herent case yet. But, on the other hand, the total vari- 
ation decays very fast and approximately remains flat, 
except for p near to 1 which in this case, as it is clear 
from figure, the total variation increases critically. It 
means that the noise applied on the position of the walker 
changes the probability distribution such that one can 
find the walker in the neighbor sites with approximately 
equal probability. But for large values of p, a gap be- 
tween probabilities of neighbor sites occurs which is very 
sensitive by applying a bit of noise. 

At the end of this section, we would like to emphasis 
that the expression in Eq. (|68|) is independent of the type 
of walking. It means that for any type of 1DQW with 
known coherent probability distribution, one can use this 
expression to find the effects of position noise on it. Sev- 
eral types of 1DQW such as biased Hadamard walk [20j . 
QW with SU(2) coin operator (50|, 1DQW with entan- 
gled coin [14|], many coins QW [5l| and so on have been 
introduced and well studied. Some aspects of specific 
types of 1DQW can be combined with the smoothness 
property of the position decoherence and opens new use- 
ful features in the study of QW. For example, Kendon et 
al. [13] have shown that the weak decoherence, both on 
the position and on the coin subspaces, can produces flat 
probability distribution where are useful in quantum in- 
formation and computation procedures. In fact, in their 
work the coin decoherence tends to change the quantum 
probability to classic one (binomial distribution around 
the origin) and the position decoherence smooths the dis- 
tribution. It is clear that applying noise on the coin sub- 
space decreases the quantumness of the system, so their 
work has made in the cost of decrease of the quantum fea- 
tures of the system such as the speed of spread and the 
variance. One of the promising interesting aspect of the 
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FIG. 7: (Color online) Total variation of probability distribu- 
tion of QW in the presence of position decoherence versus p. 
The upper panel is close up of curve for presentation of more 
details 



current work is that we could smooth the probability dis- 
tribution without the lost of the quantum features. For 
example, we can use the many coins 51| or the 1DQW 
with entangled coins [3], both having the probability 
distributions with some peaks, and smooth the proba- 
bility distribution by tuning the noise on the position 
subspace. 



IV. SUMMARY AND CONCLUSIONS 

By using the analytical expressions for the moments 
of distribution, given recently in 4J|, we have investi- 
gated the tunneling effect in the one-dimensional quan- 
tum walk. Our analytical calculation shows that the 
shifting noise (tunneling effect) on the position subspace 
of the quantum walk does not fade the quantum behav- 
ior of the system. We have brought two witnesses for our 
claim, which behave completely different when the coin 
subspace is subjected to the decoherence. First, we have 
shown that the quadratic dependency of the variance on 
the time has never vanished, even for mixed initial state. 
Second, we have calculated the entanglement between 
the coin and the position and have shown that unlike the 
coin only decoherence, in which the entanglement goes to 
zero very fast by increasing the strength of the noise, the 
entanglement in the position only decoherence (tunneling 
effect) converges to a significant value, even for maximum 
rate of noise. Furthermore, we have derived an exact ex- 
pression for the probability distribution of the QW, in 
the presence of tunneling effect, in terms of the corre- 
sponding coherent probability distributions. Our results 
have shown that the effect of the tunneling noise on the 
QW is the smoothness of the probability distribution, 
such that this smoothness occurs for all noise strength p, 
except for p near to 1 for which the smoothness is bro- 
ken. Moreover, we found that although for p near to 1, 
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the probability distributions of the even and odd sites 
are separate from each other, but even for this case the 
probabilities distribution of the odd and even sites are 
smooth individually (see Fig. [7]). Another interesting re- 
sult which we have found is that, for maximum rate of 
noise, i.e. p = 1, only the even sites are occupied despite 
of the number of steps. 

Accordingly, in the all experimental realization of the 
QW where the the quantum behavior, such as the speed 
of spreading or the entanglement, plays an important 
rule, one does not need be worry about the noise on the 
position subspace, because this type of decoherency does 
not fade the quantumness of the system, but the signif- 
icant classicality will happen whenever the coin is sub- 
jected to the decoherence. The present paper shows that 
the position only noise not only does not make the transi- 
tion from quantum to classic but also it can smooths the 
probability distribution, a task which is useful in quan- 
tum information and quantum computation process. 

Since the smoothness of the probability distribution 
and the exact expression of Eq. (TT5|) are independent 
from the type of the QW, it can be used with other as- 
pects of different type of QW in order to invent useful 



process in quantum information and quantum compu- 
tation. As an example, many coins QW introduced by 
Brun et al [5l[ , persists the quantum behavior of coherent 
1DQW for M <t, where M is the number of coins. Also 
this type of QW has many picks in the probability distri- 
bution, which can be smoothed by applying the proper 
noise on the position subspace. As an another example, 
we could mention the quantum walk on the line with 
entangled coins [14| . In this type of the QW the prob- 
ability distribution has a pick around the origin which 
can be smoothed by using the tunneling noise without 
worry about the classicality of the system. Therefore the 
results of this paper and other aspects of the lDQWs can 
be combined together and open a new window in order 
to using properties of QW in quantum information and 
quantum computation processes. 
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